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We investigate the persistence probability of a Brownian particle in a harmonic potential, which 
decays to zero at long times - leading to an unbounded motion of the Brownian particle. We con- 
sider two functional forms for the decay of the confinement, an exponential and an algebraic decay. 
Analytical calculations and numerical simulations show, that for the case of the exponential relax- 
ation, the dynamics of Brownian particle at short and long times are independent of the parameters 
of the relaxation. On the contrary, for the algebraic decay of the confinement, the dynamics at long 
times is determined by the exponent of the decay. Finally, using the two-time correlation function 
for the position of the Brownian particle, we construct the persistence probability for the Brownian 
walker in such a scenario. 



I. INTRODUCTION 

The phenomenon of persistence has been of continu- 
ing interest in the past decade. Persistence is quantified 
through the persistence probability pit), that a stochastic 
variable has not changed its sign over a time t. In a wide 
class of non-equilibrium systems this probability decays 
algebraically with an exponent and the exponent has 
been studied in systems that include free random walk in 
homogeneous [TJ E] and disordered media j3] , critical dy- 
namics [4] , surface growth [5l41 lj , polymer dynamics |12j , 
diffusive processes with random initial condition |13H15j , 
advected diffusive process [TB] and finance [T7J [TB]. A 
precise theoretical prediction for pit) can only be worked 
out only for a select few cases [19] - the simplest scenario 
being an exponentially decaying stationary correlator, as 
in the case of an overdamped Brownian motion. In gen- 
eral, for most Gaussian stochastic processes the decay of 
the stationary correlator, C(T) = (X(T)X(O)), is non- 
exponential. The behavior of C (T) in the neighborhood 
of zero characterizes the density of zero crossings for the 
underlying stochastic process [T3l [19] . When C(T) near 
zero, has a quadratic dependence on time in the first or- 
der, the number of zero crossings of the stochastic process 
is finite, and the exponent 9 is extracted using the Inde- 
pendent Interval Approximation (IIA) [13] or the sign 
time distribution of the stochastic variable [TS]. Con- 
versely, when C(T) - l-0(T a ), with a < 2, the density 
of zero crossings is infinite and perturbation expansions 
about a random walk correlator gives a good estimate of 
the persistence exponent [5J. 

The simplest of all these systems, which exhibit an al- 
gebraic decay of pit) with an exponent 1/2, is the case 
of a overdamped Brownian particle. Lying in the in- 
terface of science and engineering, Brownian motion is 
ubiquitous around us and plays a dominant role in the 
nano and mesoscopic world. The underlying principle 
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of this stochastic process is not only used for theoreti- 
cal modeling of a wide range of complex phenomena [20] , 
but Brownian motion in itself serves as an experimental 
tool for probing microscopic environments (20fl23| . In 
the popular Langevin picture, the erratic motion of a 
Brownian particle is well described by Newton's equa- 
tion of motion with a viscous drag and a delta-correlated 
stochastic force acting on the particle. While the non- 
Markovian nature of the phenomenon can be taken into 
consideration by using a generalized Langevin equation 
with a finite correlation time for the stochastic noise and 
a memory dependent friction, in the following discussion 
we shall restrict ourselves to the Markovian scenario. 

In this article, we investigate the persistence proba- 
bility of a Brownian particle in a time dependent po- 
tential - a scenario corresponding to the trapping of a 
tracer particle in some potential which eventually re- 
laxes to zero. To keep the following discussions at an 
analytically tractable level, we choose a harmonic po- 
tential, given by Uix,t) = ^fit)x 2 . The function fit) 
can be viewed as time-dependent spring constant, with 
fit) — > as t — > oo, so that the particle motion be- 
comes unbounded in the long-time limit. The converse 
situation of a constant confinement strength has already 
been studied in Ref. QH El ES] ■ In the Fokker-Planck 
description, the calculation of the persistence probability 
translates to solving the backward Fokker-Planck equa- 
tion, with an absorbing wall at the x = 0. An alterna- 
tive approach to determine the survival probability, as 
outlined in [TJ [19l El] , is from the two-time correlation 
function for the position of the stochastic variable x - ex- 
ploiting the fact that for a Gaussian stationary process 
with a correlator decaying exponentially at all times, the 
persistence probability also decays exponentially. 

The rest of the article is organized as follows: we in- 
troduce the dynamical equations of motion and construct 
the two-time correlation functions in section [ITJ A discus- 
sion on the mean-square displacement of the Brownian 
walker and the relevant time scales due to the time de- 
pendent trap is also presented in this section. We study 
two types of relaxation phenomena - an exponential and 
an algebraic relaxation of the confinement discussed in 
section |II A| and section |IIB( respectively. Finally, the 
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persistence probability is discussed in section III 



II. BROWNIAN PARTICLE IN A 
TIME-DEPENDENT POTENTIAL 

For simplicity, we take the overdamped limit for which 
the dynamics of a Brownian particle, with a unit mass, 
is governed by 



x = -f(t)x(t) + V (t), 



(1) 



where 77(f) is the stochastic velocity characterizing the 
solvent. The above equation is further supplemented by 
the moments of the stochastic noise, 

(77(f)) = and (r)(t)r){t')) = 2DS{t - t'). (2) 

At this point, we assume that the stochastic noise in "in- 
ternal" , characterized by the viscosity and the tempera- 
ture of the solvent, while the time-dependent confinement 
is "external" and does not change the delta correlation. 
An experimental realization of the model system would 
correspond to a laser trapping of a tracer, with the in- 
tensity of laser decaying in time. In such a scenario, the 
transport parameters and the temperature T gets renor- 
malized [351 ■ The corresponding solution to ([l} is 
given by 

x(t) = e-tiWW / t dt 1 77(f 1 )e/o 1 (3) 



with the initial condition x(0) = 0. Denoting g(t) — 
Jo /(Odt', the two time correlation function can be con- 
structed from ([3]), 

(x(h)x(t 2 )) = e -»(*O e -*(fe) dtj £ df 2 

Wi)n^2))e a{t ' l) e 9 ^\ (4) 

A. Exponential relaxation 

We first consider the case when the relaxation of the 
potential is given by an exponential decay, /(f) = Ae _ */ T . 
There are two time-scales in the system, r and A -1 . The 
later determines the time scale when the Brownian par- 
ticle is confined in the potential, whereas the former de- 
termines the relaxation of the potential (fig.[l£i). We fur- 
ther consider the situation when r > A -1 ; the relaxation 
time-scale of the potential is larger than the entrapment 
time-scale A -1 . When r < A -1 , the confinement decays 
even before the particle can be trapped, with the result 
that the Brownian particle undergoes free diffusion. 

Consequently, the function g(t) takes the form g(t) = 
At(1 — e -t / T ). Substituting for g(t) in Q and subse- 
quently using ([2| we arrive at 



(x(fi)x(f 2 )) = 2De 



,\ ; -e- t i /T e Are- t 2/ T 



dt' 2 e 



-2At 



(5) 



with the assumption that fi > f 2 . Performing the in- 
tegral over f 2 in ([5]), the two-time correlation function 
becomes, 



(x{t{)x(ta)) = 2Dre Xre 4l/ V Te Wt [Ei(-2A T 
- Ei(-2Are~ t2 / r )], 



(6) 



where Ei(x) is the exponential integral defined as 

/oo 
z^e-* dz. (7) 

Using ([6]), the mean-square displacement (x 2 (t)) reads 



as, 



'-(t)) = 2Dre 2XTe ~ t/T [Ei(-2Ar) - Ei(-2ATe-*/ T )] . 



. . (8) 

At this point, it is instructive to construct the limiting 
behaviors of the mean-square displacement - when t < 
A -1 and f > r. There are two scenarios we consider below 
- the first, when r is large so that the limit r — > 00 is 
appropriate, and the second when r remains finite. In the 
limit of r — > 00, the relaxation of the potential is slow 
and the Brownian particle feels a constant confinement 
strength A and (fsl) reduces to 



(x 2 (t)) 



D 



(l-e- 2 * A )+0(r- 1 ). 



(9) 



To construct the corresponding two-time correlation 
function, we expand the exponentials in ^ and keep 
the terms which are independent of t. The evaluation of 
the integral over t' 2 then gives, 



(x(ti)x(ta)) = 



,-A(ti-i a ) 



-A(ti+t 2 ) 



(10) 



which is exactly the correlation function for a non- 
stationary Ornstein-Uhlenbeck process [5S]. Eventually, 
for f >> t the particle motion becomes unbounded and 
the mean-square displacement grows linearly with time. 
A formal quantitative result in this limit can be derived 
if we take the limit of f — > 00 in Q and expand the term 
within the brackets to get 

(x 2 (f))= 2Dt e 2XTe ~ t/T [ - 7 + 2Are- 4 ^ - AVe" 24 ^ 



+ Ei(-2Ar) + - + I ln(l/2A 2 r)l . 
r 2 



(11) 



where 7 is the Euler's constant with a numerical value 
of ~ 0.5772. Keeping in mind that f >> r, the expo- 
nential functions in (11) can be ignored in comparison 



to the linearly growing term which survives, so that we 
recover the classic diffusion of the Brownian particle with 
(x 2 (t)) = 2Dt. In the opposite limit of f — > 0, a Taylor 
expansion of p| yields, 



(a; 2 (f)) = 2Dt + 0{t 2 ). 



(12) 
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The two limiting behaviors in (111 and (12) are com- 



pletely independent of the time scales, and therefore, do 
not contain any information about the confinement po- 
tential. On the contrary, when the relaxation is slow 
(t — > oo), only the short-time dynamics is independent 
of A or t. The asymptotic mean-square displacement, in 
the limit of a slow relaxation, is constant in time and is 
determined by the ratio D/X. 



B. Algebraic relaxation 

We now consider our second choice for the relaxation 
dynamics of the harmonic potential - an algebraic decay 
of the time dependent spring constant, 



(13) 



with a < 1. Using Q, the two-time correlation function 



is 



{x{h)x{t 2 )) 



2D 



du 



o/(l-o).2Aiti 



(14) 



where \\ = At q /(1 — a). The integral ( 14 ) over u yields 



2D 

(x(ti)x(t 2 )) = e 

1 — a 



(15) 



where the integral 7 is the lower incomplete gamma func- 
tion defined as 



7(0, z) = / e~ u M Q - 1 du 



(16) 



and r(a, z) is the upper incomplete Gamma function sat- 
isfying T(a, z) + 7(0, z) = T(a). The numerical value of 
7(0, z) can be evaluated using Gauss's continued fraction, 
which converges for all values of z. 

Substituting t\ = t 2 = t, the mean-square displace- 
ment is given by 

(^(t)) = J^ e -2A lt 1 -° ( _ 2Al) -l/(l- Q ) 

1 — a 



1-c* 



(17) 



Unlike the exponential relaxation of the potential, 
there is a single crossover time scale which emerges from 
(17) 



1 - a 
At" 



1/(1-") 



(18) 



and separates the regimes of normal diffusion and sub- 
diffusion in the system (fig. [Ija ) . A Taylor expansion of 
(fl7| for t < t gives, 



(x 2 (t)) = 2Dt + 0{t 2 - a ) for t<T, (19) 
while the asymptotic expansion yields 



(x 2 (t)) = 



(a^) ^ +C(^ (1_2Q) ) for t>T. (20) 



This counter intuitive result can be understood by con- 
sidering the motion of a free Brownian particle. In the 
absence of the confinement potential the Brownian par- 
ticle moves a distance y/t in time t. If we now switch 
on the potential, the strength of the potential becomes 
X(r/t) a t and for a < 1 we see that the particle feels the 
"soft" walls all the time. Mathematically, this argument 
translates to the fact that the new time scale r in ( 18 ) 



diverges as a — > 1 and is not defined in the real line for 
a > 1. We note that for a = 0, eq. (17) reduces to that 
of the Ornstein-Uhlenbeck process, 



-2Atl 



(21) 



In fig. [T] we show the mean-square displacement of a 
Brownian particle whose dynamics is governed by 
and (2}, with f(t) given by an exponential (fig. [T^,) and 
an algebraic (fig. [T]d) relaxation. The numerical integra- 
tion of ([I]) was done using the Euler scheme with an in- 
tegration time step of dt = 0.001. In the numerical solu- 
tions, the value of the diffusion coefficient D was taken as 
unity. For the exponential relaxation of the confinement, 
the measured mean-square displacement show three dis- 
tinct regimes - two diffusive regimes with a crossover in 
between. For very short (t < A -1 ) and long times (t > r), 
the particle does not feel the trap and its motion is purely 
diffusive, corresponding to (11) - (12) . In the interme- 
diate times, we observe a plateau for A -1 < t < r, corre- 
sponding to the trapping of the particle in the potential. 
To understand the origin of this plateau, we expand the 
exponential in ([T]), and retaining the zeroth order term 
then leads to a constant confinement, so that the mean- 
square displacement saturates to a value oc A -1 . This 
behavior can be observed in the left panel of fig. [TJ the 
main figure of which presents data for constant A but 
different r, while the inset shows data for a constant r 
but different A. A comparison shows that the plateau is 
determined by A -1 . On the contrary, for the algebraic 
relaxation, since a single time scale emerges from the dy- 
namics, we observe only one crossover regime determined 
by t which separates the diffusive and the sub-diffusive 
regimes (the right panel of fig.[l]). 



III. PERSISTENCE PROBABILITY 

To obtain the persistence probability, we take the route 
outlined in Ref. [TJ [TS] - we map the non-stationary 
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FIG. 1. (Color online) (a): Double logarithmic plot of the mean-square displacement (x 2 (t)} for the exponential relaxation of 
the confinement; with A = 1 and r = 1.0 (a), 10.0 (■) and 100.0 (•). The solid lines are the plots of Q for the corresponding 
values of r and A. The thick dot-dashed line in the main plot corresponds to A -1 = r = 1, while the dashed and the dot-dashed 
line corresponds to r = 10 and r = 100, respectively. Inset: We depict the variation of the mean-square displacement for 
different values of A with r = 100 fixed; A = 0.01 (□), 1.0 (O) and 10.0 (A). The solid lines are plot of |8| for the corresponding 
values of A and r. The thick dot-dashed line corresponds to A -1 = r = 100 while the dashed lines denote the values of 
A" 1 = 0.1 and A" 1 = 1. (b): Plot of mean-square displacement for A = 1.0, r = 0.001 and a = 0.1 (•), 0.2 (■) and 0.5 (A) 
and the corresponding plots of (171 for th e three values of a. The thick dashed black line is the plot of 2Dt and the thin 
dot-dashed black lines are the plots of (201 for the corresponding values of a. Inset: Plot of mean-square displacement for 



A = 1.0, a = 0.5 but with r = 0.001 (o), 0.01 (□) and 0.1 (A). 



process x(t) to a stationary Ornstein-Uhlenbeck process. 
This is usually achieved, first by a normalization of x(t) 
by \J (x 2 (t)), the root-mean-square distance the parti- 
cle has traveled and then using a suitable transforma- 
tion in time. Once, we have the stationary process X, 
with correlator C(T), the persistence problem reduces 
to calculation of no zero crossing of X. When C{T) is 
a purely exponential decay for all times, the persistence 
probability is the solution to the backward Fokkcr-Planck 
equation for an Ornstein-Uhlenbeck process, which can 
be shown to decay as P{T) = f sin" 1 [C(T)\, [U 
An application of this method, therefore, requires the 
transformation of the stochastic process x(t) to a Gaus- 
sian stationary process. Since the correlation function in 
^ and (15) is non-stationary, we make the following 



stationary Ornstein-Uhlenbeck process, 



transformations - we first define the normalized variable 
X(t) = and construct the two-time correlation 

function, following which we make a suitable transforma- 
tion in time to make the correlator stationary, as well as 
an exponentially decaying function for all times. 



Before we proceed to give a derivation of the persis- 
tence probability for the two models introduced above, 
we derive a general result applicable to the model system 
in ([!]). To transform the nonstationary process in JT]), 
we consider the transformations X = x{t)/l(t) and e — 
l 2 (t)e 2 ^\ where l 2 (t) = (x 2 (t)) and g(t) = ff(t')dt'. 
Substituting these transformations in we obtain a 



dX 



\x + rj{T), 



(22) 



where fj(T) is a Gaussian white noise with zero mean and 
unit variance. The relation between fj and r\ can be deter- 
mined from the transformation of the delta function and 
takes the form rj(T) = j§-ri(t). The stationary correlator 

for the process in ( p2| is then given by C(T) = e~ T l 2 . 
Accordingly, the persistence probability in real time de- 
cays as p(t) - e'^)/l(t). In the following, we illustrate 
this explicitly for the two specific cases presented in sec- 
tion US] and section HTBl 



A. Exponential Relaxation 

The two-time correlation function for X reads, 

(ar(*i)x(t 2 )) 



y/(x 2 (t 1 ))(x 2 (t 2 )) 

Using ([6| and ^ in the above equation, we have 



(23) 



tXU \Y(i \\ / Ei(-2Ar)-Ei(-2Are^/-) 
(X( tl )X(t 2 )) = y E . ( _ 2Ar) _ Ei( _ 2Are _ tl/T) • (24) 

The time transformation e T = l 2 {t)e 2g ^ reads, 

e T = Ei(-2Ar)-Ei(-2Are~ t/T ), (25) 
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which transforms (124]) to 



(X(T 1 )X(T 2 ))=e 



-|(Ti-T 2 ) 



(26) 



The correlator for the stochastic process X is station- 
ary and exponentially decaying. The asymptotic behav- 
ior of the persistence probability for such a process is 
then given by P(T) ~ e~ T l 2 [19 . Transforming back to 
real time, the persistence probability for the process x(t) 
is then given by, 



p(t) ~ Ei(-2Ar) - Ei(-2ATe-* /r ) 



-1/2 



(27) 



For t « A and t » r, a Taylor and an asymptotic 
expansion of the above equation gives p(t) ~ t~ x / 2 . Fi- 
nally, in the limit of r — > oo, the persistence probability 
reads as 



p(i)~ [(l- e - 2At )e 2A * + 0(r- 2 )] 



-2M-V2 



(28) 



which is identical to the result of Ref . [25 . To determine 
the persistence probability of the Brownian particle using 
a numerical integration, we chose an ensemble of random 
initial conditions in the neighborhood of zero (so that the 
sign of x(0) is well defined) and followed the sign change 
of the position. The fraction of particles which did not 
change the sign of the coordinates in time t gives an esti- 
mate of the persistence probability. The results presented 
in ( 27 ) (the colored lines in fig. [2]) and ( 28 ) (the black 
dashed line in fig. [2| is compared with the measured per- 
sistence probability using the numerical simulation of ([!]) 
(the solid points) in fig. [2j For short and long times, the 
persistence probability pit) ~ i -1 / 2 , a signature of purely 
diffusive motion presented in (11) and (12). 



B. Algebraic relaxation 

To determine the survival probability, we proceed in 
the similar way and construct the two-time correlation 
function for the normalized variable X. 



(X( tl )X(t 2 )) 



>h(t 2 ) 



(29) 



where the function h(t) is the bracketed term in (jTYJ) , 
h(t) = (-2A 1 ) 1 /( 1 -«) 7 (V(l - a), -2A 1 t 1 ^ Q ) 

Defining the time transformation e T = l 2 (t)e 29 ^ = 
h(t), the non-stationary correlator in ( 15 1 is transformed 



into a Gaussian stationary correlator which decays expo- 
nentially. Following [19] . the persistence probability, in 
real time decays as 



Pit) 



(-2A 1 )" 1 /( 1 - Q ) 7 (l/(l-«),-2A 1 t 1 ^) 



-1/2 



(30) 



a, 10 




FIG. 2. (Color online) Double logarithmic plot of persistence 
probability where f(t) decays exponentially, for A = 1.0 and 
t = 1.0 (•), 5.0 (■) and 10.0 (k). The solid and the dashed 
lines are the plots of (271 for the corresponding values of r 



and dashed line is the plot of ( 28 \ 




FIG. 3. (Color online) Plot of persistence probability for A = 
1.0, t= 0.001 a = 0.1 (•) and 0.2 (■). The solid lines in the 
main plot are plots of ( |30[ ) for the corresponding values of a. 
Inset: Plot of the persistence probability for A = 1,t = 0.01 
with a = 0.1 (•) and 0.2 (■). The solid black line is the plot 
of i -1 / 2 an d th e dashed lines are the plots of the asymptotic 
expansion ( 32 l . 



We next consider the limiting behaviors of the persis- 
tence probability given in (30). Substituting a — in 



( 30 1 , the probability reduces to the case of a harmoni- 



cally confined Brownian particle with constant confine- 
ment strength 25 , 



Pit) 



[e 2At (l 



-2At 



)] 



-1/2 



(31) 



G 



For a finite value of a < 1, when t < r, the persistence 
probability decays as p(t) ~ t" 1 / 2 while an asymptotic 
expansion of (30) gives, 



0.5 - 



Pit) 



1 



(32) 



In fig. § we compare the results of (|30|), (|31|) and (|32|) 
with the measured persistence probability from the nu- 
merical integration of (flj ) . The colored lines in the figure 
correspond to Eq. (|30[), while the dot-dashed lines are 
plots of Eq. (32). At short times, the motion is purely 
diffusive, and therefore we observe a i -1 / 2 decay of p(t) 
(the solid line in fig. [3]) . 

We note, that even though the mean-square displace- 
ment for t » t is similar to that of fractional Brownian 
motion, the decay of the persistence probabilities in the 
two scenarios are entirely different. For a particle which 
performs a fractional Brownian motion, the correspond- 
ing steady state persistence probability decays purely al- 
gebraically with an exponent 1 — a/2 [6]. 



C. Effect of Inertia 

Finally, before concluding, we remark upon the diver- 
gence of p(t) as t — > 0. This singularity is entirely the ar- 
tifact of coarse-graining in ([l ) , where we have neglected 
the inertia term. Strictly speaking, at this level of coarse- 
graining, we are not allowed to take the t — > limit, since 
the inertia of the particle plays a important role at such 
short times. The inclusion of the inertia term changes the 
short-time dynamics of the particle at to a deterministic 
one, as opposed to purely diffusive motion observed in the 
overdamped limit. Since the motion is now determinis- 
tic, the particle is persistently driven away from its initial 
position (the velocities remain strongly correlated), with 
the effect that the survival probability becomes constant. 
The purpose of this section is not only to demonstrate 
this fact that the inertia term in the Langevin equation 
indeed removes the singularity in the persistence proba- 
bility, but is also motivated by the recent experimental 
evidence of the ballistic regime of a Brownian particle 
[2"5l 121?] . While an accurate analysis would correspond to 
solving with the inertia term included, but it becomes 
difficult to extract any information from the resulting ex- 
pressions. However, since we are looking at a time much 
shorter than A -1 , exclusion of the confinement is justi- 
fied - as the particle does not feel the confinement at 
such small times. 

To this end, we consider the complete Langevin equa- 
tion for the momentum of a particle without any poten- 
tial confinement, 



P 



7 



P + r) 



(33) 



together with (rj(t)) = and (ri(t)ri(t')) = 2k B Tj6(t - 
t'). Since at short times, the dynamics of a Brownian 
particle in the model systems presented above is purely 
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FIG. 4. (Color online) The main figure depicts a double log- 
arithmic plot of persistence probability of a Brownian walker 
with the inertia term included in the Langevin equation for 
the ratio of 7/m = 0.2 (■), 1.0 (•) and 5.0 (a). The solid 
line is a plot of the approximate result derived in (371. In- 
set: Double logarithmic plot of the persistence probability 
of a Brownian particle in the overdamped limit. The solid 
line is power law fit to the data t~ e , with 6 = 0.482422. The 
analytical prediction for 9 in this limit is 1/2. 



diffusive, it suffices to consider the Langevin equation for 
a free particle for our present discussion. 

The two-time correlation function for the velocities 



decays exponentially as (v(ti)v(t2)} 
(with an initial condition (w 2 (0)) = 
sition correlation. 



ksT -\ tl -t 2 \/T 



k B T 



) and the po- 



k B T 
m 

™2 . 



+ 



2m 

1 2 - 

7 

-H/to 



-ta/ro 



-(tl-*2)/TD 



(34) 



with r = m/7 and the assumption that ti > £ 2 - While 
the transformation of the correlator in ( 34 ) to a station- 



ary process is nontrivial, we can still extract some infor- 
mation about the persistence probability in the limit of 
t — > 0. Keeping in mind this limit and that t\ > t%, a 
Taylor expansion of (34) yields, 



MtiMia)) = 



ht 2 



1*1 

2 t 



(35) 



The two-time correlation function (X(ti)X(t2)) reads as, 



(x(ti)X(t 2 ) = 



'l-*i/2r 



* 2 /2r 



(36) 



The transformation e T = (1 — t/2r ) 1 transforms (36) 
into a stationary process with an exponentially decaying 
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correlation function. The persistence probability p(t) in 
real time the translates to, 



p(t) ~ - sin" Vl - */2r ] 

7T 



(37) 



The numerical integration of ( 33 ) was done with an 



implicit integration scheme based on the Leap-Frog algo- 
rithm for different values of the ratio 7/771 and the value 
of T was set to unity. The results of the simulations 
are presented in fig. Hand compared to the approximate 
formula for p(t) in (37). 



placed by a "soft" wall. The external confinement can be 
realized using a laser-trapping experiment, with the in- 
tensity of the laser decaying in time. When the confining 
potential relaxes exponentially, we observe that the dy- 
namics of the Brownian particle at short and long times 
is purely diffusive and independent of the relaxation time 
scales. On the other hand, for an algebraic relaxation, 
the motion at long times is determined by the exponent of 
the relaxation. Using the two-time correlation function 
for the position of the Brownian particle, we construct 
the persistence probability of the Brownian particle in 
the two scenarios. 



IV. CONCLUSION 
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